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Abstract. Here we provide a detailed analysis, along with some extensions and additonal investigations, 
of a recently proposed [l] self-organised model for the evolution of complex networks. Vertices of the 
network are characterised by a fitness variable evolving through an extremal dynamics process, as in the 
Bak-Sneppen [2] model representing a prototype of Self-Organized Criticality. The network topology is 
in turn shaped by the fitness variable itself, as in the fitness network model "3" . The system self-organizes 
to a nontrivial state, characterized by a power-law decay of dynamical and topological quantities above a 
critical threshold. The interplay between topology and dynamics in the system is the key ingredient leading 
to an unexpected behaviour of these quantities. 



PACS. 89. 75. He Networks and genealogical trees - 

1 Introduction 

Complex Networks represent a very active topic in the 
field of Statistical Physics. The possibility to describe sev- 
eral different systems as structures made of vertices (sub- 
units) connected by edges (their interactions) has proved 



).65-|-b Criticality, self-organized 

successful in a variety of disciplines, ranging from com- 
puter science to biology [4, 5 . One common feature amongst 
all these different physical systems is the lack of a char- 
acteristic scale for the degree (the number of edges per 
vertex), and the presence of pairwise correlations between 
the degrees of neighboring vertices. On top of these topo- 
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logical properties, also the dynamical processes acting on These standard approaches to network modelling assume 
a network often show unexpected features. Indeed, the dy- that the topology develops much faster than other vertex- 
namics of the outbreak of an infectious disease is totally specific properties (no dynamics of the latter is consid- 
different when defined on complex networks rather than ered), or alternatively that the vertices' properties un- 
on regular lattices [6l[7] . dergo a dynamical process much faster than network evo- 
Many different models have been proposed in order to lotion (the topology is kept fixed while studying the dy- 
reproduce such behaviour. One of the earliest ones, in- namics). Of course, by separating the dynamical and topo- 
troduced by Barabasi and Albert 0, uses "growth" and logical time-scales, and allowing the "fast variables" to 
"preferential attachment" as active ingredients for the on- evolve while the "slow variables" are fixed (quenched), the 
set of degree scale invariance. On the other hand, there is ^^olc picture becomes simpler and more tractable. How- 
growing empirical evidence [HO] that networks are often ever, this unavoidably implies that the slow variables must 
shaped by some variable associated to each vertex, an as- be treated as free parameters to be arbitrarily specified 
pect captured by the 'fitness' model [8l[TTl[T2]. In general and, more importantly, that the feedback effects between 
the onset of scale invariance seems to be related to the fine topology and dynamics are neglected. 

tuning of some parameter(s) in all these models, while it The main results highlighted by the self-organised model 

would be desirable to present a mechanism through which is that these feedback effects are not just a slight modifi- 

this feature can develop in a self-organised way. nation of the decoupled scenario Rather, they play a 

. major role by driving the system to a non equilibrium sta- 
Here we describe m detail, and partly extend, a re- 

_ , • 1 tionary state with non-trivial properties, a feature shared 

cently proposed model pLj where the topological proper- 

.„ „ ... with other models of dynamics-topology coupling on net- 

ties explicitly depend on vertex-specinc fitness variables, 

, , works [HlfTillTSlfroilTT]. 

but at the same time the latter spontaneously evolve to 

self-organized values, without external fine tuning. The 

steady state of the model can be analytically solved for 2 The model 

any choice of the connection probability. For particular 

2.1 The Bak-Sneppen model and SOC 

and very reasonable choices of the latter, the network de- 
velops a scale-invariant degree distribution irrespective of order to define a self-organized mechanism for the on- 
thc initial condition. scale-invariance, we took inspiration from the activ- 

This model integrates and overcomes two separate frame- ity in the field of Self-Organized Criticality(SOC) 18J. In 

works that have been explored so far: dynamics on fixed particular, we considered a Bak-Sneppen like evolution 

networks and network formation driven by fixed variables, for the fitness variable [2 . In the original formulation. 
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one deals with a set of vertices (representing biological |24j networks. The value of the critical threshold r is found 
species) placed on a one-dimensional lattice whose edges to depend on the topology, but the stationary fitness dis- 
represent predation relationships. This system is therefore tribution always displays the same qualitative step-like 
a model for a food chain, and it has been investigated to behaviour. These more complicated structures are closer 
show how the internal dynamics is characterised by scale- to realistic food webs |25j , and allow the number of up- 
free avalanches of evolutionary events 2 . Every species i dated vertices (which equals one plus the degree of the 
is characterised by a fitness value Xi drawn from a uni- mimum-fitness vertex) to be heterogeneously distributed, 
form probability distribution p{x) defined between and Nevertheless, as long as the network is fixed, the model 
1. The species with the minimum value of the fitness is leads to the ecological paradox that, after a mutation, 
selected and removed from the (eco)system. This reflects the new species always inherits exactly all the links of 
the interpretation of the fitness as a measure of success of the previous species. This represents a problem, since it 
a particuar species against other species. Removal of one is precisely the structure of ecological connections among 
species is supposed to affect also the others. This is mod- species which is believed to be both the origin and the 
eled by removing the predator and the prey (the nearest outcome of macroevolution [26 . At odds with the model, 
neighbours in a 1-dimensional lattice) of the species with a mutated species is expected to develop different interac- 
the minimmn fitness. Three new species are then intro- tions with the rest of the system, 
duced to replace the three ones that have been removed. 

Alternatively, this event is interpreted as a mutation of 2.2 Self-organised graph formation 

the three species towards an evolved state. The fitnesses 
of the new species are always extracted from a uniform 
distribution on the unit interval. The steady state ob- 
tained by iterating this procedure is characterised by a 
uniform distribution of fitnesses between a lower critical 
threshold r = 0.66702 ± 0.00008 and 1. The size s 
of an avalanche, defined as a series of causally connected 
extinctions, follows a scale-free distribution P{s) oc s^^ 
where x = 1-073 ± 0.003 [19]. 



These problems are overcome if the network is not fixed, 
and evolves as soon as vertices mutate. At the same time, 
if the evolving topology is assumed to depend on the fit- 
ness values, one obtains a self-organised model that also 
allows to explore the interplay between dynamical and 
topological properties. 

The simplest way to put these ingredients together is 
to couple the Bak-Sneppen model P with the fitness net- 
work model [3] . In the model so obtained pj , we start by 
The behaviour of this model has also been studied on specifying a fixed number N of vertices, and by assign- 
different fixed topologies, including regular lattices pi l201 ing each vertex i a fitness value Xi initially drawn from 
121) . random graphs '22', small-world [23 and scale-free a uniform probability distribution between and 1. This 
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as rule 1: 




Fig. 1. Left: a graph at time t. The black vertex has the mini- 
mum fitness, and the two grey ones are its neighbours. Right: at 
the next time step t+1, three new fitness values are assigned to 
the vertices (light grey), and connections are established anew "n^y an amount proportional to l/fc„ where is fs 



+ =77 (1) 

where rj is uniformly distributed between and 1. In this 
case, irrespective of the number of other vertices j is con- 
nected to, Xj will be completely updated. 

Besides this rule, here we consider a weaker require- 
ment in which the fitness of each neighbour j is changed 



between the three vertices and the rest of the network. degree. We denote this choice as rule 2: 

1 kj — 1 

coincides with the Bak-Sneppcn initial state. However, kj^ ^ kj ''^^^^^ ^^-^ 

we do not fix the topology of the network, since it is de- ^^ere again 77 is a random number uniformly distributed 

termined by the fitness values themselves as the system between and 1. Rule 2 corresponds to assume that the 

evolves. More specifically, the edge between any pair of fitness xj is completely modified only if j is connected only 

vertices i and j is drawn with a fitness dependent proba- to the minimum-fitness vertex (that is, 1). If j has 

bility /(xi,2;j) as in the fitness model [3]. At each timestep f.. _ i additional neighbours, a share (fc^ - l)/kj of is 

the species with lowest fitness and all its neighbours un- unchanged, and the remaining fraction xj/kj is updated 

dergo a mutation, and their fitness values are drawn anew, to ij/kj. This makes hubs affected less than small-degree 

As soon as a vertex i mutates, the connections between it vertices. Clearly, it also implies that the probability of 

and all the other vertices j are drawn anew with the spec- connection to all other vertices varies by a smaller amount, 
ified probability f{xi,Xj). An example of this evolution 

for a simple network is shown in Fig. [T] 2.4 Arbitrary timescales for link updates 

As we now show, the stationary state of the model does 

2.3 Two rules for mutations not change even if arbitrary link updating timescales are 

introduced. In other words, the results are unchanged even 

There are various possible choices for updating the fit- if we allow each pair of vertices to be updated at any addi- 

ness of a mutating vertex. In ref. [Ij we assumed that tional sequence of tiniesteps, besides the natural updates 

each neighbour j of the minimum-fitness vertex receives occurring when any of the two vertices mutates. Two dis- 

a completely new fitness drawn from the uniform distri- tinct pairs of vertices may also be updated at different 

bution on the unit interval. We shall refer to this choice sequences of timesteps. In general, if tij denotes the time 
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when the pair of vertices i,j is chosen for update, we allow 3 Numerical results for rules 1 and 2 

a link to be drawn anew between i and j with probabil- 

In this section we present numerical results obtained by 

ity f[xi{tij), Xj{tij)]. To see that this does not change the 

simulating the model using particular choices of the mu- 

model results, we note that the fitness of each vertex i re- 

tation rule and of the connection probability. For rule 1, 

mains unchanged until it is selected for mutation, and this 

these numerical results will be confirmed by the analytical 

occurs only if i happens to be either the minimum-fitness 

results that we present later on. 

vertex or one of its neighbours. Therefore, for any pair of 

Traditionally, one of the most studied properties of the 

vertices Xi{tij) — Xi[t^j) and Xj{tij) = Xj{t^j), where 

BS model on regular lattices is the statistics of avalanches 

tij < tij denotes the time of the most recent (before tij) 

characterizing the SOC behaviour |2j. Nevertheless, as 

mutation of either i or j. Now, since at this latest update 

shown in ref. 27J, in presence of long-range [22] con- 

the links between the mutating vertex and all other ver- 

nections as those displayed by our model, the SOC state 

tices were drawn anew, it follows that t coincides with 

can be wrongly assessed in terms of the avalanche statis- 

the timestep when a connection was last attempted be- 

tics. Indeed, the absence of spatial correlations questions 

tween i and j, with probability f[xi{t^j),Xj{t[j)]. There- 

criticality even in case the avalanches are power-law dis- 

fore, between t^- and tij , a link exists between i and j with 

tributed. We then move to the study of other properties 

probability f[Xi{t^j), Xj{tij)]. On the other hand, if at time 

of the system by considering the fitness distribution p(x) 

tij the update is performed, a link will be drawn anew with 

and the degree distribution P{k) at the stationary state. 

probability f[xi[tij), Xj{tij)] = f[xi{t\j), Xj{t\j)]. There- 
fore the connection probability is unchanged by the up- 

dating event. Since the stationary state of the model only ^-l Choice of the connection probability 

depends on this probability, we find that updating events 

To to so, we first need to choose a functional form for 

do not affect the stationary state. Therefore our model is 

the connection probability f{x,y). The constant choice 

very general in this respect, and allows for rearrangements 

f{x,y) — p is the trivial case corresponding to a random 

of ecological interactions on shorter timescales than those 

graph. This choice is asymptotically equivalent to the ran- 

generated by mutations. This extremely important result 

dom neighbour variant 22J of the BS model, the average 

also means that, if the whole topology is drawn anew at 

degree of each vertex being d — p{N — 1) « pN (we drop 

each timestep, the results will be unchanged. This is a 

terms of order 1 /N from now on) . This choice is therefore 

very useful property that can be exploited to perform fast 

too simple to introduce novel effects. Indeed, the indepen- 

numerical simulations of the model. 

dence of the topology on the fitness introduces no feedback 
between structure and dynamics. 
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We therefore make the simplest nontrivial choice. If 
we require that the fitness-dependent network has no de- 
gree correlations other that those introduced by the local 
properties alone, the simplest choice for f{x,y) is p8l[29] 



f{x,y) 



zxy 



(3) 



1 + zxy 

where z is a positive parameter controlling the number 
of links. Apart for the structural correlations induced by 
the degree sequence PSII^ . higher-order properties are 
completely random, as in the configuration model [301 
[3T] . Coming back to the ecological meaning of the Bak- 
Sneppen model, the above choice is consistent with the in- 
terpretation that the more a species is connected to other 
species, the more it is fit to the whole environment: the 
larger x and y, the larger f{x,y). When z << 1, the above 
connection probability reduces to the bilinear choice 



f{x,y) = zxy 



(4) 



In this case, a sparse graph is obtained where structural 
correlations disappear. 

3.2 Stationary fitness distribution 

In Fig[2]we show the numerical results for a series of com- 
puter simulations of the model, with sizes ranging between 
iV = 100 to TV = 2000. Both rule 1 and 2 were used to up- 
date the neighbours of the minimum-fitness vertex. Some 
qualitative features can be spotted immediately. Firstly, 
all fitness values move above a positive threshold r, as in 
the traditional Bak-Sneppen model on various topologies. 
However, here the fitness distribution above the threshold 



> Rule l.N^lOO 
1 Rule 1, N^200 
: — 1 Rule 1, N ^ 500 
^ Rule 1, N^IOOO 
Rule 1, N^2000 
Rule 2, N^400 
— ■ Rule 2, N^l 




Fig. 2. Stationary fitness distribution obtained by numerical 
simulations of networks with different numbers A'' of vertices 
and using either rule 1 or rule 2. 

is not uniform. The fitness values self-organize to a scale- 
free distribution with the same exponent — 1 regardless of 
the size N of the system. This is a remarkable difference, 
and highlights the effect of the interplay between dynam- 
ics and topology. Finally, rule 1 and 2 produce a similar 
fitness distribution. In section [H where we provide an an- 
alytical solution of the model for rule 1, we demonstrate 
that the exponent is exactly —1. 

3.3 Stationary degree distribution 

In Fig. [3] we also plot the degree distribution P{k), whose 
behaviour is similar to that of p{x). In particular, we find 
that the degree distribution has a power-law shape with 
the same exponent of the fitness distribution, plus two 
cut-offs at small and large degrees. This is not surprising 
since one can derive analytically the connection existing 
between degree distribution and fitness distribution in the 
fitness model [5^155] . The expected degree k{x) of a ver- 
tex depends on its fitness x. Thus here the lower cut-off 
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Fig. 3. Stationary degree distribution obtained by numerical 
simulations of networks with different numbers N of vertices. 

corresponds to the value ^(t) and is the effect of the fit- 
ness threshold. The upper cut-off corresponds instead to 
the maximum possible value fc(l), that depends on the pa- 
rameter choice. This behaviour will be confirmed by the 
theoretical results presented in section ID 
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Fig. 4. Dependence of the critical threshold r on zN. The 
points correspond to numerical simulations (using rule 1) for 
iV = 5000 and (from left to right) z = 0.01, 0.1, 1, 10, 100, 1000. 
The solid curve is the theoretical prediction t{zN) derived in 
ea. (|17(l using the analytical solution presented in section ID 

3.5 Average fitness versus threshold 



3.4 Size dependence of the threshold 

It is clear from Figl2] that the different sizes affect the 
value of the threshold. More precisely, the latter decreases 
as system size increases. Furthermore, r also depends on 
the parameter z. To better characterise this behaviour, 
we plot in Figl4]the dependence of r on zN. We find that 
r only depends on the combination ziV, a result that we 
later confirm analytically. Indeed, the solid curve superim- 
posed to the data is obtained using the analytical results 
of section IH 



An overall measure of the evolution of the fitness values 
from the initial state to the stationary one can be ob- 
tained in terms of the average fitness (x) of vertices at 
the stationary state. In the initial state, clearly (x) — 1/2. 
As the fitness values move above the emerging thresh- 
old, (x) increases, until in the stationary state it reaches 
an asymptotic value. An important trend we identify is 
that, for fixed N, as z increases (x) decreases. Since also 
T decreases as z increases, it is interesting to monitor this 
effect by plotting (x) as a function of the threshold r. This 
is reported in figOfor z = 0.01, 0.1, 1, 10, 100, 1000. Once 
again, the simulations agree with the theoretical predic- 
tions that will be derived in the next section. 
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Fig. 5. Dependence of the average fitness (x) of vertices at 
the stationary state on r. The points correspond to numerical 
simulations (using rule 1) for A'^ — 5000 and (from right to left) 
z = 0.01, 0.1, 1, 10, 100, 1000. The solid curve is the theoretical 
prediction derived using the analytical solution presented in 
section |4l 

4 Theoretical results for rule 1 

If rule 1 is adopted, it is possible to derive complete ana- 
lytical results even for an arbitrary linking function IJ . We 
first briefly report the general analytical solution, and then 
show its particular form for specified choices of f{x,y). 

4.1 General analytic solution 

The method we used is to consider the master equation 
for the fitness distribution p{x,t) in the system at the 
stationary state [I]. If a stationary state exists then 
dp{x, t) 



dt 



for t — > cx) 



This equation can be written in the same limit as 
dp{x, t) 



dt 



r"'{x,t) -r°"\x,t) = 



(5) 



(6) 



where r*"'°"*(a;, i) denotes the fraction of vertices whose 
fitness is x entering/leaving the system respectively at 



time step t. At the steady state {t —>■ oo) these quanti- 
ties no longer depend on t, and we denote them by r'"(x) 
and r°"*(a;). It can be shown that these quantities can 
be computed separately in terms of the distribution q{m) 
of the minimum fitness. One finds 

1 i^min) 



r'"{x) 



N 



(7) 



where (kmin) = Jq q{m)k{m)dm is the expected degree 
of the vertex with minimum fitness. This relation sim- 
ply states that when the minimum is selected, then one 
vertex (the minimum itself) plus its neighbours (on av- 
erage {kmin)) are replaced by new vertices with fitnesses 
extracted from a uniform distribution. 
One can also show that 



r°"*(x) 



q{x)lN 



X < T 



P(^) Jq q{^)f{x,ni)dm X > 



(8) 



where r is defined here as the fitness value below which, in 
the large size limit, the fitness distribution p(x) = q{x)/N 
is essentially determined by the distribution of the mini- 
mum, and above which p{x) > q{x)/N, or in other words 



Npix) 

TV— »oo q[x) 



lim 



1 X < T 



(9) 



> 1 X > r 

If one now requires r"'-{x) — r°"*(a;) at the stationary 
state, it is straightforward to obtain the analytical solution 
for any form of f{x,y) [1]: 



(rA^)-i 
Pi^) = { 1 



X < T 
X > T 



(10) 



lo fi^^ m)dm 

As a remarkable result, here we find that p{x) is in general 
not uniform for x > t, in contrast with the Bak-Sneppen 
model on fitness-independent networks. 
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The value of r is determined implicitly through the ea. ([TU)l that the uniform character above r is only possi- 

normalization condition p(x)dx — 1, which implies ble when f{x,y) is independent of x and y, or in other 

dx , , , words when there is no feedback between the dynamical 

N-1 (11) 

variables and the structural properties. This highlights the 
novelty introduced by this feedback. 

The value of r depends on the topology. In particular, 
depending on how p scales with N, eq. (fTT|) leads to 



'r Jo fix, m) dm 
The topology of the network at the stationary state can 
be completely determined in terms of p(x) as in the static 
fitness model [51 [521l33j . For instance, the expected degree 
of a vertex with fitness x is k{x) = N J f{x,y)p{y)dy. 
Therefore the above solution fully characterizes both the 
dynamics and the topology at the stationary state. 1 + pN 

As in the standard BS model, in the infinite size limit 



1 



1 pN ^0 

{l + d)-^pN^d (13) 

pN ^ oo 



^ oo the distribution q{m) of the minimum fitness is These three regimes for the fitness correspond to three dif- 
uniform between and r, while almost all other values are ^rent possibilities for the topology. In particular, they are 
above r jlj In other words q{m) = 0{t — m)/T related to the percolation phase transition driven by the 

parameter p representing the density of the network. For 
4.2 Fitness-independent networks: random graphs small values of p, the graph is split into many smaU sub- 

graphs or clusters, whose size is exponentially distributed. 
The latter result generalizes what is obtained for the random- 

As p increases, these clusters become progressively larger, 

neighbour variant of the BS model [22 . Indeed, as already 
noticed, the random-neighbour model is a particular case 
of our model, obtained with the trivial choice f{x, y) — p. 
The network is a random graph independent of the fitness 
values, with Poisson-distributed degrees. We now briefly 
discuss this case as a null reference for more complicated 
choices discussed below. Our analytical result in eg. pO)) 
becomes 



and at the critical percolation threshold pc ~ 1/A^ [H 
150] the cluster size distribution is scale-invariant. When 
p > Pc & large giant cluster appears, whose size is of order 
0{N) and whose relative size tends to 1 as p ^ 1. 

Now, it is clear that below the percolation threshold 
(subcritical regime) the minimum-fitness vertex is in most 
cases isolated, with no neighbours. Thus it is the only up- 

{(rA^)^^ X < T dated vertex, and as a result all fintess values, except the 

(12) 
(prN)^^ X > T newly replaced one, tend towards 1. This explains why 

We therefore recover the standard result that, at the sta- t ^ 1 in this case. By constrast, if pN = d with finite 

tionary state, the fitness distribution is step-like, with on d > 1 (sparse regime), then a finite number of vertices is 

average one vertex below the threshold and all other val- updated, and r remains finite as A^ — > oo. This is pre- 

ues uniformly distributed between r and 1. We note from cisely the case considered in the random-neighbour vari- 
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ant |22j . that we recover correctly. Finally, if pN — > oo 
(dense regime), then an infinite number {kmin) — pN of 
fitnesses is continuously udpated. Therefore p{x) is uni- 
form between and 1 as in the initial state, and r ^ 0. 
Therefore the possible dynamical regimes are tightly re- 
lated to the topology of the network, depending on the 
parameter z. This result holds also for the more compli- 
cated case that we consider below. 

The average fitness (x) can be computed analytically 

as 



(x) = / xp{x)dx 



(14) 

Using eq. fT^ and ea. ([T^ . it is straightforward to show 
that 

1 + r 



(x) 



(15) 



so that, as expected, (x) increases from the initial value 
1/2 to a stationary value linearly dependent on t. 

4.3 Fitness-dependent networks with local properties 

We now present the analytical solution obtained with the 
simplest nontrivial choice discussed in section [5TTl f{x,y) — 
zxy/{l + zxy). With the above choice, it can be shown [T] 
that in the iV — > cx) limit eq. (|10|) becomes equivalent to 

{(riV)-i X <T 

(16) 
(tN)-^ +2/{zNt^x) x>t 

Remarkably, now p{x) is found to be the superposition 
of a uniform distribution and a power-law with exponent 
— 1. The value of t, obtained as the solution of eq. (fTT|) . 




0.001 0.005 0.010 0.050 0.100 0.500 1.000 



Fig. 6. Stationary fitness distribution p{x) for N — 5000 and 
various values of z. 



reads 



zN 



1 



zN ^0 



< ^(l){d)/d zN = d (17) 
zN oo 



As for random graphs, these three dynamical regimes are 
found to be related to an underlying topological percola- 
tion transition [1]. 

Figure [6] shows the stationary fitness distribution p{x) 
appearing in eq. (|16p . The theoretical results are in ex- 
cellent agreement with the numerical simulations shown 
previously in FiglS] Also, the theoretical dependence of r 
on zN appearing in eq. P7|) is plotted in FigH] as a solid 
line, and shown to fit the simulation results perfectly. 

We can now obtain the exact expression for the average 
fitness value {x) at the stationary state. Inserting eq.lfTB)) 
into eq. p^ we find 



{x) = -^ + ^^{1^t) = -^ + - — 5- (18) 



where in the last passage we have used the expression 
zN — — log(T^)/r^ obeyed by zN [IJ, representing the 
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inverse of eg. (fT7|) . For fixed N, eq- lfTS]) is plotted in Figl5] [^[3T] . This behaviour can be further characterized ana- 

as a function of r. Again the accordance with simulations lytically [1], by using the inverse function x{k) to obtain 

is very good. Remarkably, we find that, even if during the the degree distribution as P{k)dk — p[x{k)\dx The an- 

evolution r increases from to its asymptotic value as in alytical form of P{k) perfectly agrees with the numerical 

the case /(x, y) = p, now the stationary average value {x) results [ij. 
is not necessarily larger than the initial value 1/2. Indeed, 
for the interesting parameter range, {x) < 1/2. This is the 
effect of p{x) being no longer constant for x > t. 

5 Conclusions 



5000 




Fig. 7. Plot of k{x) for N = 5000 and various values of z. 

In this case, the expected degree k{x) — N J f{x, y)p{y)dy 
of a vertex with fitness x reads [T] 

1 -\- zx zx — ln(l + zx) 



k{x) = 



■In- 



1 + ZTX 



(19) 



This behaviour is reported in FiglTl For small (but larger 
than r) values of x, k is proportional to x. This implies 
that in this linar regime the degree distribution P{k) fol- 
lows the fitness distribution p{x), as we reported numer- 
ically in section [3731 By contrast, for large values of a; a 
saturation to a maximum degree is observed. This explains 
the upper cut-off of the degree distribution. The width of 



We have discussed in detail a recent model [T] where the 
interplay between topology and dynamics in complex net- 
works is introduced explicitly. The model is defined by 
coupling the Bak-Sneppen model of fitness evolution and 
the fitness model for network formation. The model can 
be solved analytically for any choice of the connection 
probability. Remarkably, the fitness distribution p{x) self- 
organizes spontaneously to a stationary probability den- 
sity, thus removing the need to specify an ad hoc distri- 
bution as in the static fitness model. Moreover, the sta- 
tionary state is nontrivial and differs from what is ob- 
served when dynamics and topology are decoupled. Be- 
sides providing a possible explanation for the spontanous 
emergence of complex topological properties in real net- 
works, these results indicate that adaptive webs offer a 
new framework wherein unexpected effects can be ob- 
served. 

GC acknowledges D. Donato for helpful discussions. 



the linear regime decreases as z increases, as the network This work was partly supported by the European Inte- 
becomes denser and structural correlations stronger [251 grated Project DELIS. 
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